Long Range Dynamics Related to Magnetic Impurity in the 2D Heisenberg 

Ant iferromagnet 

O.P. SushknQ 

Max- Planck- Institute fur Festkorperforschung, Heisenbergstrasse 1, D-70569 Stuttgart, Germany and 
School of Physics, University of New South Wales, Sydney 2052, Australia. 

We consider a magnetic impurity in the two-dimensional Heisenberg antifferomagnet with long 
range antiferromagnetic order. At low temperature the impurity magnetic susceptibility has a Curie 
term (oc 1/T) and a logarithmic correction (oc ln(T)). We calculate the correction and derive related 
Ward identity for the impurity-spin-wave vertex. 
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I. INTRODUCTION 

The problem of magnetic impurities interacting with a system of strongly correlated electrons has attracted much 
interest recently, mainly due to the experimental discoveries of the high-T c superconductors and new heavy fermion 
compounds. In the field of the high-T c materials, the parent compounds are known to be two-dimensional (2D) 
antiferromagnetic (AFM) Mott-Hubbard insulators based on CuC>2 planes which are driven to a superconducting 
state by doping (e.g. with holes)i£. Even though the holes can hop, thus destroying the AFM long-range order and 
causing the development of superconducting pairing, the extreme limit of static holes also has a physical relevance. 
Systems with static holes have been also realized experimentally in cuprates^ 

Several theoretical studies have addressed isolated static holes£*§iL2i£ and added spinsSiiSiii in 2D Heisenberg 
antifcrromagnets with long range AFM order. A singular logarithmic frequency behavior of the perpendicular magnetic 
susceptibility at zero temperature has been found in Ref£, see also a discussion in Ref, 12 . A very interesting problem 
is a magnetic impurity in 2D Heisenberg ant iferromagnet at 0(3) quantum critical pom^i^i^. However, this problem 
is out of the scope of the present paper. The low temperature behavior of the impurity static magnetic susceptibility 
in a gapped system is trivial, it obeys the simple Curie low, Sx = f2(f2+l)/3T, where Q is the impurity spin. However, 
for 2D systems which possess the long range AFM order at zero temperature, the excitation spectrum is gapless due 
to Goldstone spin waves and the Curie low is not obvious. A very interesting prediction^ for such a regime is the 
classical Curie low, Sx = il 2 /3T. The behaviour is classical because of alignment of the impurity moment with the 
local Neel order. This behaviour has been recently confirmed in Monte Carlo simulations for 2D S — 1/2 Heisenberg 
antiferromagnetic clusters with magnetic impurity 9 . Moreover, in these simulations a logarithmic correction, oc ln(T), 
to the classical Curie low has been found. Both the classical Curie low and the logarithmic correction are related 
to the nontrivial long-range dynamics in the system. In the present work we calculate the logarithmic correction 
using two different methods, a) Spin- wave perturbation theory, b) Semiclassical non-linear cr-model. In the leading 
1/S approximation both methods give the same result. However the results must be identical in all orders in 1/S 
and hence the comparison allows us to derive the Ward identity for the impurity-spin-wave vertex. Value of the 
logarithmic correction to the magnetic susceptibility is in agreement with Ref£. 

A crossover from quantum to classical Curie low for a finite AFM cluster with impurity is discussed in Section II. 
In Section III we derive the impurity susceptibility using the spin-wave perturbation theory, and in Section IV we 
obtain the same result using the semiclassical non-linear cr-model and derive the Ward identity. 

II. CURIE TERM 

The Hamiltonian of the system under consideration is 

H = H + H int + Hb , 
Hq = J Si ■ Sj , 

<ij> 

Hint = J±Sq ■ £1 , 



H B = B • in + ^ 



(1) 
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where Si is spin 1/2 at the site i of square lattice with antiferromagnetic interaction (J > 0), fi is the impurity spin 
coupled to the lattice spin at site 0, and B is magnetic field. To be specific we will assume that J± > 0, but all results 
are in the end independent of sign of J±. Consider anLxL cluster (L 1 is even) described by the Hamiltonian 
Hq, so no impurities for the beginning. Energy spectrum of the system is known very welU£ii£*i£. Spin of the ground 
state is 0, and lowest excitations are described by rotational spectrum of the solid top (diatomic molecule with zero 
projection of spin on axis of the molecule. K = 0) 

Es = ^±ii, (2) 

where J = 0, 1, 2... is spin of the state {J — corresponds to the ground state), / = L 2 x± is moment of inertia of 
the top, and x± ~ 0.066/J is perpendicular magnetic susceptibilityi&iii. The spectrum © is valid as soon as the 
rotation is solid, i.e. internal degrees of freedom of the top are not excited. The first internal excitation is the spin 
wave with wave length A = L (periodic boundary condition). Energy of this excitation is 

A sw = 2ttc/L, c « V2J. (3) 

There are 8 degenerate spin- wave excitations: S z = ±1, x and y directions, and two excitations (cos and sin) in each 
direction. If we consider T <C A sw then only rotations J2Jl are important. A more accurate criterion for solid rotation 
is: 8exp(-A stu /T) < 1, i.e. 

r«r—^. (4) 

In this temperature regime magnetic susceptibility of the cluster is determined by the spectrum (|2J. 

where summation over J is performed from to oo. If T <C T rot = 1/1 the susceptibility is zero, and if T rot <C T <C 
T sw the evaluation of (J5J gives 

3T \ 3 18/T y 3 A 9T V 6T y y ' 

The expansion is in integer powers of T rot /T. 

Now let us put an impurity with spin which interacts with cluster via H. intl see Eq. (Q. Excitation spectrum of 
the cluster with impurity is slightly different from Now this is a symmetric top with spin projection on axis of 
the top K = fl. This is like a diatomic molecule with uncompensated electron spin and strong spin-axis interaction. 
Rotational spectrum of such top is~ 

j(j + 1) - 2Q 2 + < n 2 > 

Ej Yl , (7) 

where J = f2, fi + 1, f2 + 2, ... is total spin of the cluster, and < i~2 2 > is the average value of J~2 2 in the intrinsic 
reference frame. Similar to the previous case the spin waves are not excited as soon as the inequality (0J is valid. 
One can easily check that in this case the magnetic susceptibility of the cluster xi is given by the same Eq. JSJ, the 
only difference is that summation over J is performed not from to oo, but from Q to oo. If T <^ T rot the cluster 
susceptibility is xi = + 1)/3T. If T rot <T< T sw evaluation of JSJ gives 

if i ~ i + 3n 2 \ 2 , i - sn 2 i + 3n 2 



Xi = — 2 IT \-fl + — + ... = -L 2 x± + — n~ + ■■■ (8) 

A 3T V 3 18/T J 3 A 9T 18/T 2 v ; 

Hence the impurity susceptibility defined as xi ~ Xo reads 

Ximp = Xi ~ Xo = ^ (l + + ■■■) ( 9 ) 

The expansion goes in integer powers of T rot /T. The leading term in 10 agrees within. 

It is interesting to note that if there are n impurities with spin f2 on the same sublattice than in Eq. © one shall 
replace — > nQ. The impurities are not independent because the cluster is rigid. In the thermodynamic limit it 
means that all the impurities within the correlation length are not independent. 
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III. LOGARITHMIC CORRECTION TO SUSCEPTIBILITY, THE SPIN- WAVE DERIVATION 

The above consideration is valid for the case when the inequality Q is valid and hence there are no real spin- wave 
excitations. However virtual spin-wave excitations are always there. In the intrinsic reference frame the cluster looks 
like picture in Fig.l 

t 

| I U If H I 

FIG. 1: Schematic picture of an AFM cluster with impurity 

The interaction Hamiltonian Hi nt , Eq. , can be rewritten in terms of spin- wave operators a and b (for review of 
the spin-wave theory see, e. g. Refi) 

H mt = J ± n-,S ^ij ± (n + 5o_ + 0-5o+) = ijx(o + 6o + 0_6j) 

q V q 

-^J^^^+M+^^+^-i^l+Vq^-q}) • (10) 

g 

Here /?| and a q are creation operators for spin waves with spin projection S z = ±1 respectively, u q and v q are 
Bogoliubov parameters, 

2 J 1 2 J 1 



and cjq '' = 2 Jyl — 7^ — * \[2Jq is the spin wave dispersion in the leading l/S approximation. In (|10f> we assume that 

S 1 = 1/2. In Eq. (|10|) the impurity-spin-wave vertex is derived in the leading 1/S approximation. However, at the 
last step we have introduced the vertex renormalization factor Zr that takes into account all higher 1/S corrections 
to the vertex. Generally speaking Zr depends on q, but here we consider only the small q limit, Zr = Z^{q = 0). 

To use perturbation theory we will assume that Ji <^ J. One loop correction to the impurity energy is given by 
diagram in Fig. 2. The diagram is infrared convergent because it contains the denominator, Ae — £| — = J± < S z >, 



FIG. 2: One loop correction to the impurity energy, the solid line shows the impurity and the dashed line shows virtual spin 
wave. 

related to the flip of the impurity spin. Here < S z >= \Z S w 0.307 is staggered magnetization of the lattice, and 
Z s w 0.61 is the renormalization factor for the staggered magnetizatio n^?'^ ^ . 
Interaction of the impurity with perpendicular magnetic held is of the form 

H B = -Bn x = -^B(n + + n_) . (12) 

The leading contribution to the impurity energy related to B is given by diagram in Fig. 3. Corresponding formula 



FIG. 3: Leading contribution to the impurity energy related to magnetic field. The dot denotes interaction with magnetic field. 
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reads 

r OS/2) 2 m , 
oeo = ■ (13) 

The contribution is finite because Ae is finite. Now let us look at one loop corrections to Seo shown in Fig. 4. These 

q 



t ; { t I i t I ; t ) ij t 



FIG. 4: One Zoop spin-wave corrections to the impurity energy related to the magnetic field. The dot denotes interaction with 
the field. 

corrections are infrared divergent because they contain the energy denominator without spin flip, ej — ep Using usual 
Schroedinger perturbation theory one finds the corresponding energy correction 



, 2 ' 
J i, 



( V2n.i ± z r \ ( V2n.j ± z r 
5e, = V I Vr^g *> + [ "> I . (14) 



Here uj q —* cq = Z c \^2Jq, where Z c w 1.17 is the spin-wave velocity renormalization due to higher 1/S 
correctionsi*i2iiS. Keeping only divergent terms we find from l|14f> the following expression for the susceptibility 

jjz 2 n 2 ^ u 2 q + v 2 q _ j\z 2 n 2 r 1 d 2 q z 2 n 2 r d q 

X± L 2 (Ae) 2 2f u q Z c (Ae) 2 J q 2 {2ir) 2 ~ 2^Z C J <S Z > 2 J q { } 

It is interesting that 114|) is independent of J±. We have to put some lower limit in the integral in (|15fl . At low 
temperature, when inequality {2} is valid, we have to substitute the spin-wave gap © as the lower limit. This gives 
temperature independent contribution to the impurity susceptibility 

7 2 2 90 2 7 2 

If J T 3> T SU) one has to put temperature as the lower limit in fTB^I . Hence 

m 2 z 2 

s ^ = ^zJj 1 ^ ™ 

In the leading 1/5 approximation one has to set Zr = Z c = Z s = 1, hence Sx±_ — > 2 |j- ln(..). This corresponds to the 
low- frequency susceptibility derived in Ref£. We have calculated the perpendicular susceptibility ill6|) . (|17fl m the 
intrinsic reference frame. The isotropic susceptibility is related to 5\± by the standard relation, S\ — § #X-L- 

IV. NON-LINEAR a-MODEL DERIVATION OF THE LOG CORRECTION, AND WARD IDENTITY 

FOR THE IMPURITY SPIN- WAVE VERTEX 

An alternative derivation of l|16(l and l|17|) is based on the a- model. Let us consider a field n, |n| = 1, defined on 
a disc of radius L. An impurity with spin f2 is in the center of the disc. The impurity spin is directed along z axis 
(perpendicular to the plane) and due to the magnetic field B = B x it is tilted by angle 9 is x-direction, see Fig. 5. 
Energy of the medium is 

E ° = \pJ (Vn)Vr , (18) 
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FIG. 5: An impurity (dot) in the center of a disc of radius L. The disc carries fild n described by the nonlinear a-model. Due 
to magnetic field B — B x the impurity spin is tilted in x-direction 

where p s = Z p j is the spin stiffness. Here i is the leading 1/S value for the stiffness, and Z p w 0.72 is the 
renormalization factor due to higher l/S'-correctionsiii2*iS. The field is of the form n « (n x ,0, 1). Due to (|I8|I the 
field n x obeys the usual Poisson equation, solution of the equation is 

n x {r)=a\n{L/r) . (19) 

To find the constant a we have to remember that n x (r *~ 1) = 9, where 9 is tilting angle of the impurity, and 
1 = lattice spacing. Therefore a = 9/\n(L). Substituting l|19|) in l|18|l we find the elastic energy E a = irp s 9 2 / ln(L). 
Total energy related to the impurity consists of the magnetic energy and the elastic energy 

Minimizing it with respect to 9 we find 9 — BQhi(L) / (2irp s ) , magnetic moment M = £19, and the magnetic suscep- 
tibility 

tt 2 20 2 

a Xx = -_ln(i)->— _ln(£). (21) 

Z7T/9 S TTZpJ 

At the final step we have substituted p s = JZ p /A. In the leading 1/S* approximation, = Z c = Z s = Z p = 1, 
hence equation l|21|) agrees with the spin- wave results (|16f) and (|17f) . Moreover a comparison of these equations gives 
a nontrivial Ward identity relating renormalization factors for the spin- wave vertex Zr, the spin- wave velocity Z c , the 
staggered magnetization Z s , and the spin stiffness Z p 

Zl - 1 (22) 



ZcZ 2 Z p 



This gives previously unknown value of Zr 



! ZcZl 



0.76 . (23) 



Similar to l|17H one has to substitute ln( J/T) instead of ln(L) in (|21|l if J ^> T > T sw . If the external magnetic field 
has a nonzero frequency u, and w > T,T SW , then ln(i) — > ln(J/o;). Equations 116() . 117|l . and l|21() are in agreement 
with Ref.- and with recent result a 9 i 21 . 

The spin-wave approach in section III assumes that Ji <^ J. There are numerous two-loop diagrams which 
are proportional to ( J±/(2tt) 2 J 2 ) ln(L) and even to ( Jj l /(2tt) 2 J 2 ) ln 2 (L). Some of the diagrams which contain the 
logarithm squared are shown in Fig. 6. Calculation of all the diagrams is a pretty involved problem. On the other 
hand the semiclassical derivation based on the a-model is independent of Jj_/J. The only assumption in the derivation 
is that the impurity magnetic moment is localized in the vicinity of the impurity. To check this assumption we have 
calculated the magnetic cloud around the impurity using the spin-wave theory. We have found that density of the 
induced magnetization Sp z drops down faster than 1/r 2 , Sp z (r) < 1/r 2 . Therefore, the magnetic moment of the cloud 
p = J Sp z (r)d 2 r is convergent at large distances, and it means that the above assumption is valid. This implies that 
all higher order in J± infrared divergent diagrams must cancel out. This is a highly unusual situation and it would 
be interesting to check the cancellation by a direct calculation. 
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FIG. 6: Two loop double logarithmic contributions to 5x± proportional to J±/J 2 . 



V. 



CONCLUSIONS 



We have analyzed excitation spectrum of a finite antiferromagnetic cluster with magnetic impurity, and considered 
a crossover between quantum and classical Curie law for the impurity magnetic susceptibility. We have also derived a 
logarithmic correction to the impurity magnetic susceptibility. Dependent on the parameters this can be logarithm of 
system size, temperature, or frequency of the external magnetic field. Using the results for the logarithmic correction 
we have derived the Ward identity for the impurity-spin-wave vertex. 
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